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Abstract 

A new dimensional scaling method for the calculation of excited states of multielectron atoms is 
introduced. By including the principle and orbital quantum numbers in the dimension parameter, 
we obtain an energy expression for excited states including high angular momentum states. The 
method is tested on He, Li, and Be. We obtain good agreement with more orthodox quantum 
mechanical treatments even in the zeroth order. 
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I. INTRODUCTION 

The dimensional scaling (D-scaling) method originally developed by Witten jlj and Her- 
schbach |2] is now a well established technique in atomic physics jjj. Substantial progress 
has been made in improving the accuracy of the technique and extending it to excited states. 
Most notable are the approaches of Goodson and Herschbach 14, nl 6f and the dimensional 

fin 

perturbation theory (DPT) developed by Carzoli, Dunn, and Watson [3, |8( . In the present 
work, an alternative approach to generalized D-scaling is developed which does not require 
calculations of high order 1/D corrections in order to obtain high-angular momentum states 
unlike, for example, the Pade approximate used by Goodson which required 20 orders of 1/D. 
Our approach is simple and yet accurate. It predicts excited states of helium (the canonical 
multi-electron problem) and is easily extended to N electron atoms. We demonstrate the 
straight forward extension by analyzing lithium and beryllium. 



II. TWO APPROACHES TO D-SCALING 

In this section, we will review the traditional D-scaling approach for helium and introduce 
our alternative quantum number dimensional formulation. 

A. Traditional D-scaling for helium 

The Hamiltonian for helium, in atomic units, is given by 

H = - 1 -(Vl + V 2 2 ) + V(r 1 ,r 2j 6), (1) 

where 

V( ri ,r 2 ,9) = -— - — + ; 1 (2) 

n r 2 ^/r\ + r\ - 2r x r 2 cos 9 

is the Coulomb potential energy, 9 is the angle between the electron radii vectors, and Z is 

the nuclear charge. The Laplacian operator, for S states, can be written in D-dimensions as 

V? + V* = Ko^n) + K d . x [t 2 ) + (\ + i J L\_ x (3) 
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where 

K o-dr) = ^T^r"'^ (4a) 

£ - = 5^»- M 4- (4b) 

m 1 2 1 2 

The probability distribution function in D- dimensions is defined to be \§d\ = Jd\^d\ 
where the D-dimensional Jacobian is given by 

Jd = (r ir2 ) D - 1 sin D - 2 6 (5) 

_ j. 

and ffl is a solution to (0) with eigenvalue E. Letting Q act on = J D 2 $ D and 
transforming parameters as 

(D-1) 2 

n -> Ti (6a) 

E - (6b) 

we obtain the D — > oo limit of the energy 

When this energy expression is minimized with respect to the parameters (ri,r 2 ,6), an 
estimate of the ground state (Is 2 ) energy of helium can be obtained which is 5% accurate. 
The accuracy can be improved by the 1/D expansion. The limitation of this approach is 
that expression (J7J) cannot produce excited S states nor higher-angular momentum states 
(states with L > 0) in its current form without resorting to Langmuir vibrations and high 
order 1/D corrections. 



B. Present D-scaling for helium 

Motivated by the approach of the previous section and to avoid high order dimensional 
perturbation expansions, we invented a generalized Laplacian of the form 

V? + V 2 , = K lx ( ri ) + K 12 (r 2 ) + ( i + 1) Ll (8) 

VI '2/ 

where 71, 72, and a are three different parameters which we also assume will tend toward 00. 
Without loss of generality, we then factor our parameters into a finite integer part multiplied 

3 



by a dimensional term given by 

71 = m(D - 1) (9a) 

72 = n 2 (D - 1) (9b) 
a = L(D-2). (9c) 

Thus, 7j and a will tend toward oo faster than D provided n, and L are finite positive inte- 
gers greater than one. By introducing the parameters n^ and L, into the usual D dimensional 
Laplacian, we will have additional degrees of freedom in the final energy expression. This 
parameterization can be interpreted as being the cardinality of a space larger than D di- 
mensions which includes the energy levels of the physical system. We will refer to this space 
as quantum number dimensional space. If we need to recover the usual three dimensional 
Laplacian, we can add a second term to Eq. given by 

7i = m(D-l) + 6(l-ni)/D (10a) 
l2 = n 2 (D-l) + Q(l-n 2 )/D (10b) 
a = L(D - 2) + 3(1 - L)/D. (10c) 

This parameterization allows us to recover the usual three dimensional Schrodinger equation 
when D = 3. But for large values of D, system El and system ITU1 are the same. In 7 — a 
space, the corresponding Jacobian factor now reads 

J 7 _ a = rfrf sin Q 6 (11) 

which is chosen to provide consistency with the traditional approach. If we now solve 
H^rf-a, = Ety^^ a with definitions (|5jl. (|9"|). and the scaling transformation © we obtain the 
final energy expression in the D — > 00 as 

The derivation of this expression is given in the Appendix A. Additionally, in the Ap- 
pendix B we show how to correct the energy with the first 1/D correction. Equation (|12|) will 
allow us to obtain atomic excited states by associating nj as the principle quantum number 
for electron i and L as the sum of the individual hydrogenic assignment orbital quantum 
numbers i.e. L = t\ + 1 2 . Thus, our L is the maximum in the set of possible total orbital 
angular momenta L 6 {\£ x — £ 2 \, . . . , \£± +-^2|}- 
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TABLE I: Energy (in atomic units) of the He atom obtained using Eq. (|13j) for stats with L = 0. 
The final column is the percent error between our energy and previously published results. 



\n\n 2 L) 



Assignment 



-E (Present.) 



-E (Ref.) a 



A(%) 



|110) 
1 120) 
1 130) 
1 140) 



Is 2 !S 
ls2s 1 S 
ls3s 1 S 
IsAs l S 



3.0625 
2.1595 
2.0621 
2.0333 



2.9037 
2.1460 
2.0613 
2.0336 



5.40 
0.60 
0.03 
0.02 



"Complex rotation technique 
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III. GROUND STATE AND EXCITED S STATES 



The excited S states of helium can easily be found by setting L = in expression (fT2"|) . 
we find 

E = l(ji + ^)+V(n,r 2 ,e). (13) 

It is interesting to note that this energy expression is identical to what one would obtain 
for a Bohr model of helium using the Bohr space quantization condition for the kinetic 
energy operators i.e. Inr = nX where A is the de Broglie wavelength A = h/p. Bohr space 
quantization was also used by Greenspan [^] to obtain an energy expression for the excited 
states of helium but with a potential function that was parametrized. More recently, the 
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Here, we take the 



Bohr model has been successfully applied to molecules 
electron-nuclear and electron-electron interaction to have their usual Coulombic form. The 
results for excited S states are presented in Table ID We see that a dominant feature is the 
increase in accuracy with larger quantum number n 2 . This can easily be explained in terms 
of Bohr correspondence. Physically, as n 2 > rii{— 1) the corresponding radii will be very 
different since rj ~ n 2 . Thus, electron two sees a screened nucleus of charge Z = 1 and 
therefore it is a hydrogenic problem, which is exactly soluble in the Bohr picture. 

IV. EXCITED P AND D STATES 



Our goal here is to stay within the geometric configuration picture originally established 
by Hershbach and co-workers. Other groups have also made substantial contributions (see 
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TABLE II: Energy of the He atom for states with L > 0. 



\n\n 2 L) 



Assignment 



-E (Present.) 



-E (Ref.) 



1 121) 
1 131) 
1 132) 
1 221) 
1 222) 
1231) 



ls2p 1 P 
ls3p 1 P 
ls3d 1 D 
2s2p 1 P 
2p 2 1 D 
2s3p 1 P 



2.12710 
2.05573 
2.05572 
0.69720 
0.68440 
0.55990 



2.12380 & 
2.05514 & 
2.05562 fe 
0.69310 c 
0.68641 c 
0.56384 c 



0.20 
0.020 
0.004 
0.60 
0.29 
0.70 



b Reference 
c Reference 



m 
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TABLE III: Positions of the electrons on the hypersphere for certain representative states 



\nin 2 L) 



n 



1 121) 
1 132) 
1 222) 
1231) 



0.5005 
0.5000 
2.4278 
2.0249 



3.8981 
8.9579 
2.4278 
7.9728 



90.4526° 
90.0221° 
95.3006° 
95.8483° 



[^1 and references therein) to use D-scaling for these states. We present here, for the first 
time, a true Lewis structure picture for helium excited states thus providing an example of a 
classical analogue for P states. Physically, the P and D states are just special configurations 
of the two electron system in hyperspherical coordinates. Table ITT1 compares the energy of 
the He atom for some states with L > obtained in the present analysis with the accurate 
values known in the literature. Table 11111 shows the configurations of the electrons for some 
representative states. We have found that pure S states admit an angle of 180° whereas states 
with P character and higher- angular momenta typically lie in the range of 90° < 9 < 180°. 
We cannot, with this method, distinguish between singlet and triplet states of helium nor 
can we determine the correct spectroscopic assignment of L. We note that in general, our 
method finds a value somewhere near the average of the various L states. For example, the 
2^2 ige energy is ~ —0.62 a.u. whereas for the 2p 2 3 P e state E pa —0.71 a.u. and 
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the average of the two yields —0.67 a.u. which is much closer to our result of —0.68 a.u. 
However, our calculated energy is only 0.29% different from the 1 D state. This suggests 
that our L is indeed the sum of the individual t± s. 

Our advantage over the dimensional perturbation treatment is twofold. Firstly, in order 
to calculate the excited states of helium they needed to invoke harmonic oscillator wave- 
functions as an ansatz to the problem. The various excited modes of the normal coordinates 
corresponded to the hydrogenic excited states of helium. Their harmonic oscillator quantum 
numbers where associated with the radial and orbital nodes of the hydrogenic wavefunctions 
(see for example I n or der to calculate with this method, many orders of 1/D n (in the 

1/.D expansion) where required. Furthermore, it was seen that the 1/D series was asymptotic 
and required further re-summation via Pade approximates. The second drawback was the 
fact that it could only consider states where l\ = £2- Thus, the ls2p state is not accounted 
for in their treatment. 

The states with electron 1 in an s-orbit and electron 2 in a higher angular momentum 
state seem to work the best. We can correct the errors by considering small vibrations 
about the "frozen" configuration positions. The derivation of the first 1/D correction and 
its application are provided in the Appendix B for the excited states of helium. 

V. GENERALIZATION TO THE N-ELECTRON SYSTEM 

The present approach admits a simple generalization to the N electron atom with nuclear 
charge Z. We can describe the N electron system by including all pairwise interactions. 
The problem will thus involve (1/2)N(N + 1) parameters. Let 




(14) 



as before and 




(15) 



then 




(16) 



N . N 



Here 7, is as before and a^- = Lij(D — 2) with Ljj = £j + £j. Next we transform the N 
particle wave-function as 

N . N 

I*) = J[r~^ JJsi7T^%|$). (17) 

Transforming in this way makes the analysis a natural continuation of the helium prob- 
lem with the exception that the angular momentum are multiply counted. Then in the 
large-D limit the energy function to be minimized is 

1 N n 2 

E{n 1 i 1 n 2 l2 ■ ■ ■ n N £ N ) = - -|+ 

4 (iv^i) E 4 + ^) cot2 % + vtf *5, (is) 

where 

^(r, 6) = -Z ~ + + r ? _ 2riT i cos ^i)"'- ( 19 ) 

* i<j 

The factor of l/(iV — 1) accounts for the multiple counting of each li. For the case of the 
ground state of lithium [Z = 3, N = 3) in the state ls 2 2s we obtain E = —7.7468 which 
is 3.5% different from the true ground state of —7.4780. However, the accuracy of higher 
states (including li > 0) is typically much better (see table II VJ) . Additionally, we tested the 
method on beryllium and had some success. In this study, we considered states with the core 
electrons fixed in the Is 2 orbits. We found an error with respect to published values of 2.6 
percent. To eliminate this systematic error, we subtracted off the ground state error from 
each of the excited states and found very good agreement with published results. Shown 
in the table are some calculated excited states of beryllium which are compared to other 
published results. 

VI. CONCLUSION 

We have introduced an alternative form of dimensional scaling for the excited states of 
atoms. By treating the principle and orbital quantum numbers to be innately coupled to the 
dimension D, we developed an algebraic equation for multiple excited states of multielectron 
atoms including states of high orbital momentum. This method can be extended to include 
relativistic effects for nuclei with high Z numbers by including the kinetic-energy mass 
correction term of the Breit-Pauli Hamltonian. 
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TABLE IV: Results for various excited states of Lithium 





Assignment 


-E (Present.) 


-E (Ref.) 


A(%) 


|101021) 


ls 2 2p 2 P 


7.6897 


7.4101 d 


3.70 


|102021) 


Is2s2p 4 P 


5.2715 


5.4375 e 


3.10 


|102121) 


ls2p 2 l D 


5.1982 


5.2335 / 


0.70 


(102130) 


Is2p3s X P 


5.0744 


5.0820^ 


0.15 


(102131) 


Is2p3p 1 P 


5.0668 


5.0651-f 


0.03 


|102141) 


ls2pAp 3 S 


5.0395 


5.0311 / 


0.16 


(202032) 


2s 2 3d 3 P 


1.9496 


1.9690* 


1.00 


(212132) 


2p 2 3d 3 P 


1.8192 


1.8246 ft 


0.30 



^Reference 
e Reference 
■^Reference 
^Reference 
''Reference 



IS 

m 
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TABLE V: Ground state and various excited states of beryllium. Core electrons Is 2 are implicit 
in the notation. The percent error is wrt the corrected value column. 





Assignment 


-E(Present.) 


Corrected 


-E(Ref.)* 


A(%) 


(2020) 


2s 2 l S 


15.0603 


14.6673 


14.6673 


0.00 


(2021) 


2s2p l P 


14.8939 


14.5009 


14.4734 


0.20 


(2030) 


2s3s 1 S 


14.8014 


14.4084 


14.4182 


0.07 


|2121) 


2p 2 l D 


14.7674 


14.3744 


14.4079 


0.23 


(2031) 


2s3p l P 


14.7693 


14.3763 


14.3931 


0.11 


(2032) 


2s3d X D 


14.7682 


14.3752 


14.3737 


0.01 
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APPENDIX A: SCHRODINGER EQUATION IN "D DIMENSIONS" 



Here we demonstrate the idea of our D-scaling transformation by applying it to the He 
atom. The Hamiltonian of He in three dimensions is given by Eq. (JIJ. First we perform a 
continuous transformation of the Laplacian as follows 



V 2 + V 2 



1 d 



„7i 



d 



+ 



1 d 



„1'2 



d 



rj 1 dri 1 dri ^ rj 2 dr 2 2 dr 2 V r i r \) sin a 9 d9 



Id d 

(Al) 



89 



where 71, 72 and a are given by Eqs. 0. The wave function, coordinates and the energy 
transform according to 

(D-l) 2 * 



m = [rfrf sin a #n 1/2 $ 



The radial part of the Laplacian acting on \I/ yields 



E 



(D-l) 



:E. 



1 d 



r 7 Or \ dr J 
while the angular part leads to 

sm a 9d9\ 89 



= r 2 s in 2 x 



dr 2 



a 

39 2 + 2 



7^ ii 
2, 



a\ 1 a 

2. 



= r 2 sin 2 ^ x 

After the transformation the Schrodinger equation reduces to 

(T + U + V)& = E<$>, 



cot 2 9 



where 



T 



2 


- d 2 d 2 




32 


(£> - l) 2 


dr\ dr 2 


yr 2 + r 2 J 


d9 2 



is the kinetic energy term and 

_ 7i(7i-2) 1 _ 72(72 ~ 2) 1 _ / 1_ ^ _l_ 



2(D-l) 2 r 2 2(D-1) 2 



a (a — 2) o 



r 2 ; L2P-1) 2 



core 



a 



(D-iy 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



is the centrifugal potential. In the large-D limit the kinetic energy term vanishes and Eq. 
()A5|) reduces to finding the minimum of the algebraic energy function (j!2j) . 
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APPENDIX B: CALCULATION OF 1/D CORRECTION 



Let us take into account terms of the order of 1/D in Eq. (|A5|i . then it reduces to 



2 

D 2 



1 



"? r 2 . / <9(9 2 



a 2 



1 /n 2 n 2 \ L 2 ( 1 1 . , t . 



i 



^ + 5+(l + l)(L(L + l)cot^ + L)|U = ^ (Bl) 

L'l '2 V'l '2/ J J 

We decompose the effective potential {7 + V in Eq. (|B1J) near the minimum (rio, T2o ; #o) an d 
leave only terms quadratic in displacement from this point. As a result, we obtain 





1 

D 


h + 


r 20 


u + 


r 2 J 






2 


Fd 2 


<9rf 








D 2 




[rio 



+ 



20 



de 2 



[AiAr 2 + /? 22 Ar 2 + (3 33 A6 2 + 2/3i 2 AnAr 2 + 2/3i 3 AnA0 + 2/3 23 Ar 2 A#]$ = (B2) 

where are coefficients of the Taylor expansion of the function U + V at the minimum. 
Eq. (jB2|) describes three coupled one dimension harmonic oscillators. To make the oscillator 
masses equal we rescale AO as AO = A0/r o , where l/r = a/1/V 2 + l/r 2 . Then Eq. (jB2[) 
yields 



. ' 10 ' 20 V ' 10 ' 20 / J J ^ 



r d 2 



+ 



d 2 



+ 



d 



2 1 



dr\ dr\ dO 2 



[/5nAr 2 +/3 22 Ar 2 2 +/333A^Vro 2 +2/3 12 Ar 1 Ar 2 +2/3 1 3Ar 1 A0/ro+2/3 23 Ar 2 A^/r o ]^ = E<$>. (B3) 

The problem is reduced to determining eigenvalues Ai, A 2 , A3 of the symmetric matrix: 

f Pn P12 /WO 

P12 P22 /W r o 
^/W r o /W r o (3 33 /r 2 j 



E 



(D-iy 



1 


m 


D 






D . 



+ + 3") (i(i + Dcot 2 «o + i) 

.'10 '20 V10 '20/ 



+ 



(B4) 
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TABLE VI: Improvement after the first 1/D correction. 



State 


Zeroth order 


First order 


Is 2 !S 


5.40% 


1.67% 


2s3p X P 


0.70% 


0.54% 



We have applied Eq. (jB4|) to correct the energy we found for the Is 2 state from table and 
the 2s3p state from [HJ Shown in table IVII is a comparison between the zeroth order and 
first 1 /D correction for these states with respect to the more exact values. We have noticed 
that it is possible to correct the energy with the 1 /D expansion to first order but only if the 
zeroth order error is significant. For states where the zeroth order gave very good accuracy 
(ie j 1%) we do not expect much improvement from the first 1/D correction. It was noted 
in reference j^] that 30 terms of the 1/D expansion leads to a 9 decimal place accuracy for 
the ground state energy of helium. Thus, the first order term should only correct the most 
significant figure. 
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